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Abstract 



In this series of lectures, presented at the CIMPA Winter School on Discrete 
Integrable Systems in February 2003, we give a review of the application of Lie 
point symmetries, and their generalizations, to the study of difference equa- 
tions. The overall theme could be called "continuous symmetries of discrete 
equations" . 
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1 Introduction 



1.1 Symmetries of Differential Equations 

Before studying the symmetries of difference equations, let us very briefly review 
the theory of the symmetries of differential equations. For all details, proofs 
and further information we refer to the many excellent books on the subject 
e.g. |lHl[71|inil31,22,3,51inH|. 

Let us consider a completely general system of differential equations 

Ea{x,U,Ux,Ua:x, ■ ■ -Unx) ^ 0, X G , U G , a = 1 , . . . iV, (1) 

where e.g. Unx denotes all (partial) derivatives of order n. The numbers p, q, n 
and N are all nonnegative integers. 

We are interested in the symmetry group G of the system (^, i.e. in the 
local Lie group of local point transformations taking solutions of eq. into so- 
lutions. Point transformations in the space X xU of independent and dependent 
variables have the form 

x^Ax{x,u), u = Qxix,u), (2) 
where A denotes the group parameters. We have 

Ao{x,u) — X, flQ{x,u)—u 

and the inverse transformation {x,u) — > exists, at least locally. 

The transformations 10) of local coordinates in X x U also determine the 
transformations of functions u — f{x) and of derivatives of functions. A group 
G of local point transformations oi X x U will be a symmetry group of the 
system ^ if the fact that u{x) is a solution implies that u{x) is also a solution. 

The two fundamental questions to ask are: 

1. How to find the maximal symmetry group G for a given system of equa- 
tions (UJ? 

2. Once the group G is found, what do we do with it? 

Let us first discuss the question of motivation. The symmetry group G allows 
us to do the following. 

1. Generate new solutions from known ones. Sometimes trivial solutions can 
be boosted into interesting ones. 

2. Identify equations with isomorphic symmetry groups. Such equations may 
be transformable into each other. Sometimes nonlinear equations can be 
transformed into linear ones. 

3. Perform symmetry reduction: reduce the number of variables in a PDE 
and obtain particular solutions, satisfying particular boundary conditions: 
group invariant solutions. For ODEs of order n, we can reduce the order 
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of the equation. In this reduction, there is no loss of information. If we 
can reduce the order to zero, we obtain a general solution depending on 
n constants, or a general integral (an algebraic equation depending on n 
constants) . 

How does one find the symmetry group G? One looks for infinitesimal trans- 
formations, i.e. one looks for the Lie algebra L that corresponds to G. Instead 
of looking for "global" transformations as in eq. ^ one looks for infinitesimal 
ones. A one-parameter group of infinitesimal point transformations will have 
the form 

x.i = Xi + X^i{x,u) |A| << 1 (3) 
Ua = Ua + \(j)a {x,u) I < i < p, ^ < a < q. 

The functions and (jja must be found from the condition that u{x) is a solution 
whenever u{x) is one. The derivatives Ua,xi must be calculated using eq. lj2J) 
and will involve derivatives of and 4'a- A K-th derivative of Ua with respect 
to the variable Xi will involve derivatives of £,i and 4>a up to order K. We then 
substitute the transformed quantities into eq. and request that the equation 
be satisfied for u{x), whenever it is satisfied for u{x). Thus, terms of order A° 
will drop out. Terms of order A will provide a system of determining equations 
for and Terms of order A'^, fc = 2, 3, . . . are to be ignored, since we are 
looking for infinitesimal symmetries. 

The functions and (j)a depend only on x and u, not on first, or higher 
derivatives, Ua,Xi, Ua.XiX^j ^tc. This is actually the definition of "point" sym- 
metries. The determining equations will explicitly involve derivatives of Ua, up 
to the order n (the order of the studied equation) . The coefficients of all linearly 
independent expressions in the derivatives must vanish separately. This provides 
a system of determining equations for the functions £,i{x,u) and 4>a{x,u). This 
is a system of linear partial differential equations of order n. The determining 
equations are linear, even if the original system is nonlinear. This "lineariza- 
tion" is due to the fact that all terms of order A^ , j > 2, are ignored. 

The system of determining equations is usually overdetermined, i.e. there are 
usually more determining equations than unknown functions and (jja (p + Q 
functions). The independent variables in the determining equations are x G R^, 
u G R«. 

For an overdetermined system, three possibilities occur. 

1. The only solution is the trivial one = 0, = 0, « = 1, . . .p, a = 1, . . . , g. 
In this case the symmetry algebra is L — {0}, the symmetry group is G = J 
and the symmetry method is to no avail. 

2. The general solution of the determining equations depends on a finite 
number K of constants. In this case the studied system has a finite- 
dimensional Lie point symmetry group and we have dim G = K . 

3. The general solution depends on a finite number of arbitrary functions of 
some of the variables {x^, Uq}. In this case the symmetry group is infinite 
dimensional. This last case is of particular interest. 



3 



The search for the symmetry algebra L of a system of differential equations is 
best formulated in terms of vector fields acting on the space XxU oi independent 
and dependent variables. Indeed, consider the vector field 

p 9 
X = '^(_i{x,u)dxi + ^^(l)a{x,u)dua, (4) 

1=1 Q=l 

where the coefficients and are the same as in eq. ■ If these functions are 
known, the vector field Q can be integrated to obtain the finite transformations 
©. Indeed, aU we have to do is integrate the equations 

^=^,(x,'u), ^ = 0a(x,M), (5) 
subject to the initial conditions 

Xi\x=o=Xi Ua\\=O^Ua. (6) 

This provides us with a one-parameter group of local Lie point transformations 
of the form Q with A the group parameter. 

The vector field Q tells us how the variables a; and u transform. We also 
need to know how derivatives like u^, Uxx, ■ ■ ■ transform. This is given by the 
prolongation of the vector field X. 

We have 

prX^X + Y,{T. ^'o:du., + E ^T'du^.x, (7) 

;,fc 



XiXkXl N 

i,k,l ' 



where the coefficients in the prolongation can be calculated recursively, using 
the total derivative operator 

^Xi ^Xi ^~ ^a.Xi^Uct ~^ ^a.XaXi^^Ua-Xa ~^ ,X aXf^X i -X aX ^ "t" . . . (8) 

(a summation over repeated indices is to be understood). 
The recursive formulas are 



4>a' = Dj:^(f)a - {Dx^S,a)Ua,x^ 
r^;^" = D^,r^^ - {D^,^a)Uo.,x.x. (9) 
(bT-"" = D^.ct^T" - {D,,^aW,,^x,x. 

etc. 

The n-th prolongation of X acts on functions of x, u and all derivatives of 
u up to order n. It also tells us how derivatives transform. Thus, to obtain 
the transformed quantities we must integrate eq. (|SJ with conditions ©, 
together with 

= Uj), Ui\x=0=Ux. (10) 
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We see that the coefficients of the prolonged vector field are expressed in terms 
of derivatives of and (/iq,, the coefficients of the original vector field. They carry 
no new information: the transformation of derivatives is completely determined, 
once the transformations of functions are known. 

The invariance condition for the system is expressed in terms of the 
operator as 

pr(")Xi;jB,^._£„^o=0, a=l,...N, (11) 

where pr*^"-' X is the prolongation Q calculated up to order n where n is the 
order of the system 

In practice the symmetry algorithm consists of several steps, most of which 
can be carried out on a computer. For early computer programs calculating 
symmetry algebras, see Ref. [5^1101 • For a more recent review, see |25) . 

The individual steps are: 

1. Calculate all the coefficients in the n-th prolongation of X. This depends 
only on the order of the system i.e. n, and on the number of indepen- 
dent and dependent variables, i.e. p and q. 

2. Consider the system ([3) as a system of algebraic equations for x, u, Ux, 
Uxx, etc. Choose N variables wi, V2, ■ ■ - Vn and solve the system ^ for 
these variables. The Vi must satisfy the following conditions. 

(i) Each Vi is a derivative of one Uq, of at least order 1. 

(ii) The variables Vi are all independent, none of them is a derivative of 
any other one. 

iii) No derivatives of any of the Vi figure in the system (|T|l. 

3. Apply pr^"^ X to all the equations in ^ and eliminate all expressions Vi 
from the result. This provides us with the system 1)11(1 . 

4. Determine all linearly independent expressions in the derivatives remain- 
ing in Hll|) . once the quantities Vi are eliminated. Set the coefficients of 
these expressions equal to zero. This provides us with the determining 
equations, a system of linear partial differential equations of order n for 
4>a(x,u) and^i(a;,u). 

5. Solve the determining equations to obtain the symmetry algebra. 

6. Integrate the obtained vector fields to obtain the one-parameter subgroups 
of the symmetry group. Compose them appropriately to obtain the con- 
nected component Go of the symmetry group G. 

7. Extend the connected component Go to the full group G by adding all 
discrete transformations leaving the system invariant. These discrete 
transformations will form a finite, or discrete group Gd- We have 

G = Gd ^ Go (12) 

i.e. Go is an invariant subgroup of G. 
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Let us consider the case when at Step 5 we obtain a finite dimensional Lie 
algebra L, i.e. a vector field X depending on K parameters, K £ Z>, K < oo. 
We can then choose a basis 

{XuX2,...,Xk} (13) 

of the Lie algebra L. The basis that is naturally obtained in this manner depends 
on our integration procedure, though the algebra L itself does not. It is useful 
to transform the basis (|13() to a canonical form in which all basis independent 
properties of L are manifest. Thus, if L can be decomposed into a direct sum 
of indecomposable components, 

Li(SL2®...®Lm, (14) 

then a basis should be chosen that respects this decomposition. The components 
Li that are simple should be identified according to the Cartan classification 
(over C) or the Gantmakher classification (over R) |47II24| . The components 
that are solvable should be so organized that their nilradical |52l32j is manifest. 
For those components that are neither simple, nor solvable, the basis should be 
chosen so as to respect the Levi decomposition |52ll32j . 

So far we have considered only point transformations, as in eq. lj2Jl, in which 
the new variables x and it depend only on the old ones, x and u. More general 
transformations are "contact transformations" , where x and u also depend on 
first derivatives of u. A still more general class of transformations are generalized 
transformations, also called "Lie-Backlund" transformations |48,^. For these 
we have 

X Ax{x,U,U:c,Ua;x, ■ ■ ■) (15) 

U = nx{x,U,Ua:,Uxx,---) 

involving derivatives up to an arbitrary order. The coefficients and (jfa of the 
vector fields Q will then also depend on derivatives of Ua ■ 

When studying generalized symmetries, and sometimes also point symme- 
tries, it is convenient to use a different formalism, namely that of evolutionary 
vector fields. 

Let us first consider the case of Lie point symmetries, i.e. vector fields of the 
form Q and their prolongations ((Tj). To each vector field I0J we can associate 
its evolutionary counterpart X^ , defined as 

Xe = Qa{x,U,Ux)dUa, (16) 
Qa = K-i.^. (17) 

dxi 

The prolongation of the evolutionary vector field is defined as 

prXe = Qadua + Ql'dua.x, + QT^dua.x.x^ + • ■ ■ (18) 
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The functions Qa are called the characteristics of the vector field. Notice that 
Xe and prXe do not act on the independent variables Xi. 

For Lie point symmetries evolutionary and ordinary vector fields are entirely 
equivalent and it is easy to pass from one to the other. Indeed, eq. H17|) gives 
the connection between the two. 

The symmetry algorithm for calculating the symmetry algebra L in terms 
of evolutionary vector fields is also equivalent. Eq. Hll() is simply replaced by 

pr(")Xe£;a |Bi=...=i;„=o,= 0, a^l,...N. (19) 

The reason that eq. Hll|) and (|19|) are equivalent is the following. It is easy to 
check that we have 

pr(") X, = pr(") X - i,D,. (20) 

The total derivative Di is itself a generalized symmetry of eq. (0, i.e. we 
have 

A^^a |Bi=B.=...=is„=o,= i p, a=\,...N. (21) 

Eq. (EOJ) and prove that the systems ifTT)) and lfTO|l are equivalent. Eq. if^ 
itself follows from the fact that DEa = is a differential consequence of eq. (^), 
hence every solution of eq. 1^ is also a solution of eq. (|^ . 

To find generalized symmetries of order k we use eq. (|16ll but allow the char- 
acteristics Qa to depend on all derivatives of up to order k. The prolongation 
is calculated using eq. H18|l . The symmetry algorithm is again eq. (|19|l . 

A very useful property of evolutionary symmetries is that they provide com- 
patible flows. This means that the system of equations 

is compatible with the system In particular, group invariant solutions, i.e. 
solutions invariant under a subgroup of G are obtained as fixed points 

Qa = 0. (23) 

If Qa is the characteristic of a point transformation then H23|) is a system of 
quasilinear first order partial differential equations. They can be solved, the 
solution substituted into and this provides the invariant solutions explicitly. 



1.2 Comments on Symmetries of Difference Equations 

The study of symmetries of difference equations is much more recent than that 
of differential equations. Early work in this direction is due to Maeda |44II45 | 
who mainly studied transformations acting on the dependent variables only. A 
more recent series of papers was devoted to Lie point symmetries of differential- 
difference equations on fixed regular lattices |^|lTl|^|gg||gl|ggl|lgl|5niinilgl ■ 
A different approach was developed mainly for linear or linearizable difference 
equations and involved transformations acting on more than one point of the 
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lattice |2(JII21ll^ll28ll36| . The symmetries considered in this approach are really 
generalized ones, however they reduce to point ones in the continuous limit. 

A more general class of generalized symmetries has also been investigated 
for difference equations, and differential-difference equations on fixed regular 

lattices uniiniii^Esi- 

A different approach to symmetries of discrete equations was originally sug- 
gested by V. Dorodnitsyn and collaborators [ni[n|[niEl|Tni|Tl[Tn|[Tll|SllIZ|- 

The main aim of this series of papers is to discretize differential equations while 
preserving their Lie point symmetries. 

Symmetries of ordinary and partial difference schemes on lattices that are 
a priori given, but are allowed to transform under point transformations, were 
studied in Ref. |33|2H1|1S1. 



2 Ordinary Difference Schemes and Their Point 
Symmetries 

2.1 Ordinary Difference Schemes 

An ordinary differential equation (ODE) of order n is a relation involving one 
independent variable x, one dependent variable u — u{x) and n derivatives u, 
u, ... 

i?(x,u,H,u,...,M(")) = ^^0. (24) 

An ordinary difference scheme (OAS) involves two objects, a difference equa- 
tion and a lattice. We shall specify an OAS by a system of two equations, both 
involving two continuous variables x and u{x), evaluated at a discrete set of 
points {xn}- 

Thus, a difference scheme of order K will have the form 

Eai{x,y^:!^j,,AukKt:!+M) - 0, a = 1, 2 

K = N - M + 1, n, M, N eZ, Uk = u{xk). 

At this stage we are not imposing any boundary conditions, so the reference 
point Xn can be arbitrarily shifted to the left, or to the right. The order K of 
the system is the number of points involved in the scheme H25|) and it is assumed 
to be finite. We also assume that if the values of Xk and Uk are specified in 
[N — M) neighbouring point, we can calculate their values in the point to the 
right, or to the left of the given set, using equations H25|l . 

A continuous limit for the spacings between all neighbouring points going to 
zero, if it exists, will take one of the equations H25I) into a differential equation 
of order K' < K, another into an identity (like = 0). 

When taking the continuous limit it is convenient to introduce different quan- 
tities, namely differences between neighbouring points and discrete derivatives 
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like 



U ■ 



1_L 'y "JAi 



(26) 



-^n — 1 

In the continuous limit, we have 



^0, h_ —t 0, u^j; — s-w, u_,j. — s-u, M , 



As a clarifying example of the meaning of the difference scheme ()25|l. let 
us consider a three point scheme that will approximate a second order linear 
difference equation: 



Un +l " 2ttn + Un - 1 
E2 = Xn+1 - 2x„ + Xn-1 = 0. (28) 



^1 7Z TTi " ^' (-^^) 



The solution of eq. E2 — Q, determines a uniform lattice 

Xn = hn + xq (29) 

The scale h and the origin xq in eq. 129|l are not fixed by eq. H28|l . instead they 
appear as integration constants, i.e. they are a priori arbitrary. Once they are 
chosen, eq. (1271 reduces to a linear difference equation with constant coefhcients, 
since we have x„+i — Xn — h. Thus, a solution of eq. (|27|l will have the form 

u„-A"". (30) 

Substituting into lf?7|l we obtain the general solution of the difference 
scheme (I2Z|), (EHJ as 



u{Xn) = CiA^"+C2A2", Xn = hn + Xo, (31) 



The solution H31|) of the system H27|l . H28|l depends on 4 arbitrary constants ci, 
C2, ft. and Xq. 

Now let us consider a general three point scheme of the form 

Ea{Xn-l,Xn,Xn+l:Un-l,Un:Un+l) ^ 0, a=l,2 (32) 

satisfying 

detf^^^l^O, det 0, (33) 



d{Xn+l,Un+l) J \d{Xn-l,Un-l) 
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(possibly after an up or down shifting). The two conditions on the Jacobians 
H33II are sufBcient to ahow us to calculate (a;„+i,u„+i) if (a;„_i, m„_i, a;„, m„) 
are given. Similarly, (a;„_i, it„_i) can be calculated if (a;„, u„, a;„+i, m„+i) are 
given. The general solution of the scheme H32|) will hence depend on 4 arbitrary 
constants and will have the form 

Un = /(a;n,Cl,C2,C3,C4) (34) 
Xn = (/)(n, Ci,C2,C3,C4). (35) 

A more standard approach to difference equations would be to consider a 
fixed equally spaced lattice e.g. with spacing h = 1. We can then identify the 
continuous variable x, sampled at discrete points Xn, with the discrete variable 
n: 

Xn = n. (36) 
Instead of a difference scheme we then have a difference equation 

mukY^^^,,) = 0, (37) 

involving K — N — M + 1 points. Its general solution has the form 

Un = f[n,ci,C2,. ■ -Cn-m) (38) 

i.e. it depends on — M constants. 

Below, when studying point symmetries of discrete equations we will see the 
advantage of considering difference systems like the system (|25|l . 



2.2 Point Symmetries of Ordinary Difference Schemes 

In this section we shall follow rather closely the article . We shall define the 
symmetry group of an ordinary difference scheme in the same manner as for 
ODEs. That is, a group of continuous local point transformations of the form 
(jJJ taking solutions of the OAS H25|l into solutions of the same scheme. The 
transformations considered are continuous, and we will adopt an infinitesimal 
approach, as in eq. We drop the labels i and a, since we are considering 
the case of one independent and one dependent variable only. 

As in the case of differential equations, our basic tool will be vector fields of 
the form In the case of OAS they will have the form 

X = ^[x,u)d^+c^{x,u)du (39) 

with 

X = Xni U = Un — ^(•^n)' 

Because we are considering point transformation, ^ and cf) in l|39|l depend on x 
and u at one point only. 
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The prolongation of the vector field X is diff'erent than in the case of ODEs. 
Instead of prolonging to derivatives, we prolong to all points of the lattice fig- 
uring in the scheme ()25|l . Thus we put 

n+N n+N 

prX^ ^ ^ixk,Uk)dx^+ ^ (j)[xk,Uk)du^. (40) 

k=n+M k=n+M 

In these terms the requirement that the transformed function u{x) should satisfy 
the same OAS as the original u{x) is expressed by the requirement 

wXEa\E^=E^=o=Q. a =1,2. (41) 

Since we must respect both the difference equation and the lattice, we have 
two conditions from which to determine ^(x,u) and (j){x,u). Since each 
of these functions depends on a single point (x, u) and the prolongation H40|l 
introduces N — M + 1 points in space X x U, the equation H41|l will imply a 
system of determining equations for ^ and (j). Moreover, in general this will be 
an overdetermined system of linear functional equations that we transform into 
an overdetermined system of linear differential equations mi2| ■ 

To illustrate the method and the role of the choice of the lattice, let us start 
from a simple example. The example will be that of difference equations that 
approximate the ODE 

(42) 

on several different lattices. 

First of all, let us find the Lie point symmetry group of the ODE ()42|l . i.e. 
the equation of a free particle on a line. Following the algorithm of Chapter ^ 
we put 

pr(2) X = + (pdu + Vd, + </)""a, (43) 

U U 

= D,^(j) - {D^C) U= + (0« - ^x) U -iuU^ 
(j}^"" = D^(t>=' - {D^O 4'xx + - ^xx) u 

+ {(l)uu ~ '2.ixu)u^ - iuuU^X + (0„ - 2^^) U 

— 3^„w uu . 

The symmetry formula (|ll|l in this case reduces to 

|„^^= 0. (44) 

Setting the coefficients oi u^ , m^, u and (u'^) equal to zero, we obtain an 8 
dimensional Lie algebra, isomorphic to sl(3,R) with basis 

Xi = dx, X2^xdx, Xs^udx (45) 

X4 = du, X5 = xdu, Xe^udu, 

X7 = x{xdx + udu), X^ = u{xdx + udu)- 
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This result was of course already known to Sophus Lie. Moreover, any 
second order ODE that is linear, or linearizable by a point transformation has a 
symmetry algebra isomorphic to sl(3, R). The group SL(3, R) acts as the group 
of projective transformations of the Euclidean space E2 (with coordinates x, u)). 

Now let us consider some difference schemes that have eq. H42|l as their 
continuous limit. We shall take the equation to be 

- 2u„ + Un^i _ ^ ^^g^ 

However before looking for the symmetry algebra, we multiply out the denom- 
inator and investigate the equivalent equation 

El = Un+l - 2un + u„-i = 0. (47) 

To this equation we must add a second equation, specifying the lattice. We 
consider three different examples at first glance quite similar, but leading to 
different symmetry algebras. 

Example 1. Free particle (|47|l on a fixed uniform lattice. We take 

E2 = Xn ~ hn - xq = 0, (48) 

where h and Xq are fixed constants (that are not transformed by the group (e.g. 
h = 1, Xq = 0). 

Applying the prolonged vector field (|in|) to eq. we obtain 

C(x„,u„) = (49) 

for all Xn and w„. Next, let us apply (|40(l to eq. 147() and replace x„, using H48|l 
and Un+i, using l(T7|) . We obtain 

(f>{h{n + 1) + Xq, 2Un — Un-l) — 2(j)(hn + Xq, Un) 

+<j){h{n-l) + xo,Un-i) ^0. (50) 

Differentiating eq. H50() twice, once with respect to u„_i, once with respect to 
Un, we obtain 

0(x„+i,u„+i) = (51) 

OUn+l 

and hence 

)=A{xn)un + B{xn). (52) 

We substitute eq. (|52|l back into 1)50(1 and equate coefficients of w„, w„_i and 1. 
The result is 

A{n + l)^A{n), B{n + 1) - 2B{n) + B{n - 1) ^ 0. (53) 
Hence we have 

A^Ao, B = Bin + Bq = bix + bo (54) 
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where Aq, Bi, Bq, bi and bg are constants. We obtain the symmetry algebra of 
the OAS H47|l . (|48|l and it is only three-dimensional, spanned by 

Xi = du, X2 = xdu, X3 = udu- (55) 

The corresponding one parameter transformation groups are obtained by inte- 
grating these vector fields (see eq. 

Gi : X ~ X 

u{x) ~ u{x) + A 
G2 ■ X ^ X (56) 

u{x) = u{x) + Xx 
G3 : X — X 

u{x) — e^u{x) 

Gi and G2 just tell us that we can add an arbitrary solution of the scheme to 
any given solution, G3 corresponds to scale invariance of eq. (|47|l . 

Example 2. Free particle (I47|l on a uniform two point lattice. 
Instead of eq. H48|) we define a lattice by putting 

E2 = Xn+l - Xn = h, (57) 

where /i is a fixed (non-transforming) constant. Note that H57|l tells us the 
distance between any two neighbouring points but does not fix an origin (as 
opposed to eq. I0H|))- 

Applying the prolonged vector field (|4()|l to eq. H57(l and using (I57|l , we obtain 

^{Xn + h,Un+l) - ^{Xn,Un) ^ 0. (58) 

Since Un+i and ii„ are independent, eq. (|58|l implies ^ ^ ^{x). Moreover ^(a;„ -I- 
h) = £_{x) so that we have 

^ = ^0 = const . (59) 

Further, we apply pv X to eq. ()47|l . and put = 2w„ — x^+i — x^ + h, 

Xn-i ~ Xn — h in the obtained expressions. As in Example^we find that ^(x, u) 
is linear in u as in (|52|l and ultimately satisfies 

(/)(x, u) = au + bx + c. (60) 

The symmetry algebra in this case is four-dimensional. To the basis elements 
H55|) we add translational invariance 

Xi^d,. (61) 

Example 3. Free particle (I47II on a uniform three-point lattice. 
Let us choose the lattice equation to be 

E2 = Xn+l - 2Xn + Xn-1 = 0. (62) 
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Applying pr X to E2 and substituting for x„+i and Un+i, we find 

^{2xn - Xn-i,2un ~ - 2^(a;„,M„) +^(a;„_i,u„_i) = 0. (63) 

Differentiating twice with respect to u„ and Un-i, we obtain that ^ is linear in 
u. Substituting ^ — A{x)u + B{x) into Ht)3|) we obtain 



where A, . . ., are constants. Finally, we obtain a six-dimensional symmetry 
algebra for the OAS (|T7|l . (jHlJl with basis Xi, . . ., as in eq. It has been 

shown that the entire sl(3,R) algebra cannot be recovered on any 3 point 



^From the above examples we can draw the following conclusions. 

1. The Lie point symmetry group of an OAS depends crucially on both 
equations in the system H25|) . In particular, if we choose a fixed lattice, as 
in eq. (|48|l (a "one-point lattice") we are left with point transformations 
that act on the dependent variable only. 

If we wish to preserve anything like the power of symmetry analysis for 
differential equations, we must either go beyond point symmetries to gen- 
eralized ones, or use lattices that are also transformed and that are adapted 
to the symmetries we consider. 

2. The method for calculating symmetries of OAS is reasonable straightfor- 
ward. It will however involve solving functional equations. 

The method can be summed up as follows 

1. Solve equations H25|) for two of the quantities entering there, to make the 
equations explicit. For instance, take the system We can solve 

e.g. for Xn+i and Un+i and obtain 



2. Apply the prolonged vector field to eq. (^5)) and substitute for 
Xn+l^ Un+i- We obtain two functional equations for ^ and (j) of the form 



£,{Xn,Un) = AUn + Bx,, + C. 

Similarly, applying pr AT to eq. l|T7|l . we obtain 

(t>{Xn,Un) = DUn + EXn + F. 



(64) 



(65) 



OAS. 



(66) 




= a = 1,2 



dX', 



n— 1 



'a 



(67) 
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3. Assume that the functions (j>, Ei and E2 are sufficiently smooth and 
differentiate eq. (|67() with respect to the variables Xk and Uk so as to obtain 
differential equations for ^ and (f). If the original equations are polynomial 
in all quantities we can thus obtain single term differential equations form 
(|67|l . These we must solve, then substitute back into (|67|l and solve this 
equation. 

We will illustrate the procedure on several examples in Section 

2.3 Examples of Symmetry Algebras of OAS 

Example 4. Monomial nonlinearity on a uniform lattice. 
Let us first consider the nonlinear ODE 

u-it^ = 0, N^0,1. (68) 

For N 7^ —3 eq. ()68|l is invariant under a two-dimensional Lie group, the Lie 
algebra of which is given by 

= 9„ X2 = iN- l)xd^ - 2udu (69) 

(translations and dilations). For N = —3 the symmetry algebra is three- 
dimensional, isomorphic to sl(3,R), i.e. it contains a third element, additional 
to ()59fl . A convenient basis for the symmetry algebra of the equation 

u = (70) 

is 

Xi= dx, ^2 = 2xdx + udu, X3 = x{xdx + udu). (71) 
A very natural OAS that has (|68|l as its continuous limit is 

^ + „^A-^0 ^_^o,l (72) 

\Xn-]-l Xjij 

E2 = Xn+1 ~ 2Xn + Xn-1 = (73) 

Let us now apply the symmetry algorithm described in Chapter 12.21 to the 
system (|72() and (|73|l . To illustrate the method, we shall present all calculations 
in detail. 

First of all, we choose two variables that will be substituted in eq. (|41|) . once 
the prolonged vector field l|in|l is applied to the system ifT^ and (|75|l , namely 

Xn+i = 2a;„ - Xn-i (74) 

Un+l = (Xn - X„„i)^U^ -|- 2u„ - 

We apply pr X of (gDJ) to eq. (O and obtain 

^(a;„+i, u„+i) - 2^(.T„, Un) + f (x„_i, u„_i) = (75) 
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where, a;„, u„ Xn-i, Wn-i are independent, but Un+i are expressed in 

terms of these quantities, as in eq. 1741) . Taking this into acccount, we differen- 
tiate (|73|) first with respect to m„-i, then Un- We obtain successively 

{N{Xn- Xn-lfu^'^ +i\tur,+ iUr, + i{Xn+l,Un+l) = 0. (77) 

Eq. (|77|l is the desired one-term equation, ft imphes 

£,{x,u) = a{x)u + b{x) (78) 

Substituting (0 into (O we obtain 

- a(2a;„ - Xn-i) -t- a(a;„_i) = 0. (79) 

Differentiating with respect to a;„, we obtain a = oq = const. Finally, we 
substitute H78|) with a = ao into H75|l and obtain 

a = 0, 6(2x„-a;„_i)-25(a;„)-f 5(a;„_i) = (80) 

and hence 

^^b^bix + bo (81) 

where 6o and hi are constants. To obtain the fimction (/)(x„,u„), we apply prX 
to eq. (|72|l and obtain 

i) - 2(/)(a;„,u„) -f 0( 

-{Xn - Xn^lf[N(f>{Xn,Un)u^-^ + 2blU^] - 0. (82) 

Differentiating successively with respect to and u„ (taking (|74ll into ac- 

count) we obtain 

- 4>,u„ + iiXn+i,Un+l) + 4>,u„-i{Xn-l,Un-l) = (83) 

{7V(a;„-a;„_i)2u^ + 2}0,„„+,„„+, = (84) 

and hence 

if) = if)iu + (f)Q(x), 01 = const . (85) 

Eq. (|82|l now reduces to 

0o(2x„ - x„_i) - 2(/)o(a;„) -f (?!)o(a;„_i) 

-(x„ - x„_i)2{(Ar _ 1)01 + 2&i}u^ 

-iV(x„ - xn-ifc^ou^-^ = 0. (86) 

We have N 0,1 and hence H86|l implies 

00 = 0, (iV - 1)01 + 26i = 0. (87) 

We have thus proven that the symmetry algebra of the OAS (|72|l and H73|l is 
the same as that of the ODE (jSHl, namely the algebra 
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We mention that the value iV = — 3 is not distinguished here and the system 
1)7211 and H73|l is not invariant under SL(3, R) for N = —3. Actually, a differ- 
ence scheme invariant under SL(3,R) does exist and it will have eq. (|70|l as its 
continuous limit. It will however not have the form (|71|) and the lattice will not 
be uniform |16II18| . 

Had we taken a two-point lattice, Xn-^-i — Xn — h with h fixed, instead of 
-E2 = as in H73|l , we would only have obtained translational invariance for the 
equation (|72|1 and lost the dilational invariance represented by X2 of eq. (|69|) . 

Example 5. A nonlinear OAS on a miiform lattice 

^ un+r-2u^ + u ^_^u^-u^.,\^^^ ^^^^ 



E2 = Xn+1 - 2Xn + Xn-1 = 0, (89) 

where f{z) is some sufficiently smooth function satisfying 

/ (z) ^ 0. (90) 
The continuous limit of eq. H88|l and (|89|l is 

u -/(u) = (91) 
and is invariant under a two-dimensional group with Lie algebra 

^1=5,, X2 = du (92) 

for any function f{u). For certain functions / the symmetry group is three- 
dimensional, where the additional basis element of the Lie algebra is 

— {ax + bu)dx + {cx + du)du- (93) 

The matrix 

^-i: \ ) (94) 



c 



can be transformed to Jordan canonical form and a different function f{z) is 
obtained for each canonical form. 

Now let us consider the discrete system H88|l and ()89|l . Before applying pr X 
to this system we choose two variables to substitute in eq. (|41f) . namely 

Xn+l = 2x„ - Xn-1 (95) 

Un+1 = 2Un - Un-1 + (Xn ~ Xn-l) J [ 

\Xn Xfi — 1 

Applying pr X to eq. H89|l we obtain eq. H75|l with Xn+i and as in eq. (|95|) . 

Differentiating twice, with respect to and u„ respectively, we obtain 

^,«„+in„ + i [1 + {Xn - 2:„_l)/'][2 + (X„ - X^-l)/'] + tu„ + J" = 0. (96) 
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For /" 7^ the only solution is C.m„+i = 0; i-C- C = (,{^)- Substituting back into 
(f75|) . we obtain 

^ = ax + (3 (97) 

with a ~ const, /3 = const. 

Now let us apply pr X to Ei of eq. and and replace x„+i, Un+i as 
in eq (|95|l . We obtain the equation 

0(x„+i,u„+i) - 2(/)(x„,u„) + 0(a;„_i,u„_i) = 2a(x„ - a;„_i)^/(z) 

+ (x„ - .„_i)^ /'(^) f _ az) (98) 

with a as in eq. H97|) . Thus, we only need to distinguish between a = and 
a = \. Eq. H98|) is a functional equation, involving two unknown functions <j) 
and /. There are only four independent variables involved, a;„,a::„_i,u„ and 
Un-i- We simplify H98|l by introducing new variables {x,u,h,z}, putting 

Xn = X, Xn+i = X + h, Xn-1 = X — h (99) 
Un = U, Un-l—U—hz, Un+1 = U + hz + f (z) , 

where we have used eq. (|95|l and defined 

Un ~ Un-1 , , , 

Z= h = Xn+l-Xn- (100) 

•^n ^n — 1 

Eq. H98|l in these variables is 

(j)(^x + h,u + hz + f{z)) — 2(j)(x, u) + 0(a; — h,u — hz) 



2ah'^f{z) + h^f'iz) 



{x, u) — (j){x — h,u — hz) 



az 



(101) 



First of all, we notice that for any function f(z) we have two obvious sym- 
metry elements, namely Xi and X2 of eq. (|92|) . corresponding to a — 0, P — 1 
in (|101|l (and ^7^) and </) = and 0=1, respectively. Eq. (|101|l is quite difficult 
to solve directly. However, any three-dimensional Lie algebra of vector fields in 
2 variables, containing {Xi,X2} of eq. as a subalgebra, must have X3 of 
eq. (|93|) as its third element. Moreover, eq. (|97|l shows that we have & = in 
eq. (ESIl and In (|101|l we put a = a and 

0(x, m) = cx -|- c?w. (102) 

Substituting into eq. (|101|l we obtain 

{d - 2a)f{z) = [c+{d- a)z]f'{z). (103) 

^From eq. (|1U3|) we obtain two types of solutions: 
For d ^ a we have 

/ = /o[(rf-a)z-l-c](''-2-)/(^~'^), c^O. (104) 



18 



For d — a, we have 

/ = /oe-('^/^)^ (105) 

With no loss of generaUty we could have taken the matrix (|94|) with & = 
to Jordan cannonical form and we would have obtained two different cases, 
simplifying (|104|) and (|105|l . respectively. They are 

iV — 2 

/ = /o^^, = xd, + j;^udu, N^l (106) 

/ = /oe-^ X:i ^ xd.^ + {x + u)du. (107) 

The result can be stated as follows. The OAS (|88() and (|89|l is always invariant 
under the group generated by {Xi, X2} as in It is invariant under a three- 
dimensional group with algebra including X3 as in eq. if / satisfies eq. H103|l . 
i.e. has the form (|106ll . or H1U7|) . These two cases also exist in the continuous 
limit. However, one more case exists in the continuous limit, namely 

U= (1 + (^)2)3/2gfearctan; (^Qg) 

with 

X:i ^ {kx + u)d.:, + [ku - x)du. (109) 

This equation can also be discretized in a symmetry preserving way (Ifi) . not 
however on the uniform lattice H89|) . 



3 Lie Point Symmetries of Partial Difference Schemes 
3.1 Partial Difference Schemes 

In this chapter we generalize the results of ChapterElto the case of two discretely 
varying independent variables. We follow the ideas and notation of Ref. |SH1- 
The generalization to n variables is immediate, though cumbersome. Thus, 
we will consider a Partial Difference Scheme (PAS), involving one continuous 
function of two continuous variables u[x,t). The variables (a;,t) are sampled 
on a two-dimensional lattice, itself defined by a system of compatible relations 
between points. Thus, a lattice will be an a priori infinite set of points Pi lying 
in the real plane R^. The points will be labelled by two discrete subscripts Pm,n 
with — oo<m<oo,— oo<n<oo. The cartesian coordinates of the point Pmn 
will be denoted {xmmtmn) [or similarly any other coordinates {amm Pmn)]- 

A two-variable PAS will be a set of five relations between the quantities 
{x, t, u\ at a finite number of points. We choose a reference point Pmn = P 
and two families of curves intersecting at the points of the lattice. The labels 
m = rriQ and n = uq will parametrize these curves (see Fig^. To define an 
orientation of the curves, we specify 

•^m+ln •^m,n — ^m ^ 0, ^mn — ^ (-^-^0) 

at the original reference point. 
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X 



Figure 1: 



The actual curves and the entire PAS are specified by the 5 relations 

1 < a < 5 ii <i <i2 ji < j < 32- 

In the continuous limit, if one exists, all five equations Hlll(l are supposed 
to reduce to a single PDE, e.g. Ei — can reduce to the PDE and Ea = 0, 
a > 2 to = 0. The orthogonal uniform lattice of Fig. [3 is clearly a special case 
of that on Fig. [T] 

Some independence conditions must be imposed on the system e.g. 



1^1 



9(^1,..., 



d(^Xm-\-i2 jfi 1 t'm-\-i2ni Xm,n--\-j2 i ^m,Ti+j2 i ^m+Z2 ,^,4-^2 ) 



^0. (112) 



This condition allows us to move upward and to the right along the curves 
passing through Pm,n- Moreover, compatibility of the five equations Hlll|) must 
be assured. 

As an example of a PAS let us consider the linear heat equation on a uniform 
and orthogonal lattice. The heat equation in the continuous case is 

ut^u^x- (113) 

An approximation on a uniform orthogonal lattice is given by the five equa- 
tions 

= h2 ^ ^ 
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m + I 



Figure 2: 



E2 — Xm+l^n — Xm,7i — /ll — £-3 — im+l,n ~ tmji — (US) 

E4 — — Xmn =0 = tm,n+l — tm,n — /l2 = 0. (116) 

Equations H115|l can of course be integrated to give the standard expressions 

Xmn^him + XQ Unn = h2n + tQ. (117) 

Notice that hi and /12 are constants that cannot be scaled (they are fixed in 
eq. (|115|l . On the other hand (xo,to) a-re integration constants and are thus 
not fixed by the system l|115(l . pi6|l . As written, these equations are invariant 
under translations, but not under dilations. 

Finally, we remark that the usual fixed lattice condition is obtained from 
pi6|) by putting xq = to = 0, hi = h2 ~ 1 and identifying 

X = m, t = n. (118) 

Though the above example is essentially trivial, it brings out several points. 

1. Four equations are indeed needed to specify a two-dimensional lattice and 
to allow us to move along the coordinate lines. 

2. In order to solve the PAS (|114|1 . H115|) for hi and /i2 given, we must 
specify for instance {xmn, tmn,Umn, Um+i,n, Um-in}- T hen w e can directly 
calculate {xm+i,n,tm+in}, {xmn+1 , tm,n+i} from eq. H115|). In order to 
calculate the coordinates of the fourth point figuring in eq. H114|) , namely 
{xmn-ii tm.n~i} wc must shift cq. (I115|l down by one unit in m. 

3. The Jacobian condition (|112|l allowing us to perform these calculations, is 
obviously satisfied, since we have 

d{Ei, E2, E2, E4, E5) 



^(^m+l,n ; ^m+l,n ; '^m.n-fl : ; ''^mn+l ) 



= 1. 



(119) 
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A partial difference scheme with one dependent and n independent variables 
will involve + 1 relations between the variables (xi,a;2, . . .a;„,w), evaluated 
at a finite number of points. 

3.2 Symmetries of Partial Difference Schemes 

As in the case of OAS treated in Chapter [2 we shall restrict ourselves to point 
transformations 

X = F\{x,t,u) t — G\{x,t,u), u = H\{x,t,u). (120) 

The requirement is that ux(x, t) should be a solution, whenever it is defined and 
whenever u{x,t) is a solution. The group action (|f 20|l should be defined and 
invertible, at least locally, in some neighbourhood of the reference point Pmn, 
including all points Pm+i,n+j involved in the system Ullll) . 

As in the case of a single independent variable we shall consider infinitesi- 
mal transformations that allow us to use Lie algebraic techniques. Instead of 
transformations H12U|) we consider 

X = X + \£_{x,t,u), 

t = t + \T{x,t,u), (121) 
u = u + X(f){x,t,u) |A| << 1. 

Once the functions ^, r and (p are determined from the invariance require- 
ment, then the actual transformations H12U|) are determined by integration, as 
in eq. |(SJ), ®. 

The transformations act on the entire space {x,t,u), at least locally. This 
means that the same function F, G and H in eq. (|120|l . or ^, r and (j) in eq. (|121|l 
determine the transformations of all points. 

We formulate the problem of determining the symmetries H121(l . and ulti- 
mately (|120|) . in terms of a Lie algebra of vector fields of the form 

X — S^{x, t, u)dx -I- t{x, t, u)dt + 4'{x, t, u)du, (122) 

where ^, t and 4> are the same as in eq. (|121|) . The operator (|122|l acts at one 
point only, namely {x,t,u) = {xmn,tmn,Umn)- Its prolongation will act at all 
points figuring in the svstem l|lll|l and we put 

prX = Y^[^{xjk,tjk,Ujk)d^^^ + T{xjk,tjk,Ujk)dt^^ (123) 

+4'{xjk,tjk,Ujk)dujJ, 

where the sum is over all points figuring in eq. (|lll|l . To simplify notation we 
put 

S.jk = ^{Xjk,tjk,Ujk), Tjk = T{Xjk,tjk,Ujk) (124) 
(t)jk = (f>{Xjk,tjk,Ujk)- 
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The functions ^, r, and </) figuring in cq. (|122ll and p23|l are determined 
from the invariance condition 



pr XEa \ei=...=E5=o^ a=l,...5. 



(125) 



It is eq. (|125ll that provides an algorithm for determining the symmetry 
algebra, i.e. the coefficients ^, t and cj). 

The procedure is the same as in the case of ordinary difference schemes, 
described in Chapter |21 In the case of the system Hlll|) . we proceed as follows: 

1. Choose 5 variables Va to eliminate from the condition H125|) and express 
them in terms of the other variables, using the system and the 

Jacobian condition (|112|l . For instance, we can choose 



The quanties Va must be chosen consistently. None of them can be a 
shifted value of another one (in the same direction). No relations between 
the quantities Va should follow from the system Once eliminated 

from eq. (|124|l . they should not reappear due to shifts. For instance, the 
choice (|126f) is consistent if m + 12 and n + j2 are the highest values of 
these labels that figure in eq. Hlll|) . 

2. Once the quantities Va are eliminated from the system (|125|) . using H127() . 
each remaining value of Xi^k, ii,fe and Ui^k is independent. Each of them 
can figure in the corresponding functions ^ik, nt, (j)ik (see eq. (|124ll '). in 
the functions Ea directly, or via the expressions Va, in the functions ^, 
T and (j) with different labels. This provides a system of five functional 
equations for ^, r and (j). 

3. Assume that the dependence of ^, r and (jj on x, t and u is analytic. 
Convert the obtained functional equations into differential equations by 
differentiating with respect to xtk , Uk , or uik ■ This provides an overdeter- 
mined system of differential equations that we must solve. If possible, use 
multiple differentiations to obtain single term differential equations that 
are easy to solve. 

4. Substitute the solution of the differential equations back into the origi- 
nal functional equations and solve these. The differential equations are 
consequences of the functional ones and will hence have more solutions. 
The functional equations will provide further restrictions on the constants 
and arbitrary functions obtained when integrating the differential conse- 
quences. 

Let us now consider examples on different lattices. 




(126) 



(127) 
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3.3 The Discrete Heat Equation 



3.3.1 The Continuous Heat Equation 

The symmetry group of the continuous heat equation (|113|l is well known |48| . 
Its symmetry algebra has the structure of a semidirect sum 

L = Lo55Li, (128) 

where Lq is six-dimensional and Li is an infinite dimensional ideal representing 
the linear superposition principle (present for any linear PDE). A convenient 
basis for this algebra is given by the vector fields 

Po = dt, D ^ itdt + 2xd^ + udu 

K = 'it{tdt + xd^) + {x^ + 2t)udu (129) 

Pi = dx,B = 2tdx + xudu, W = udu 

S = S{x,t)du, St-S^^^O. (130) 

The sl(2,R) subalgebra {Pq, D, K} represents time translations, dilations and 
"expansions". The Heisenberg subalgebra {Pi,B,W} represents space transla- 
tions, Galilei boosts and the possibility of multiplying a solution u by a constant. 
The presence of 5* sjrnply tells us that we can add a solution to any given solu- 
tion. Thus, S and W refiect linearity, Pq and Pi the fact that the equation is 
autonomous (constant coefficients). 



3.3.2 Discrete Heat Equation on Fixed Rectangular Lattice 

Let us consider the discrete heat equation l|114l) on the four point uniform 
orthogonal lattice (|115|l . (|116l) . We apply the prolonged operator H113|) to the 
equations for the lattice and obtain 

^(^m,n+l 7 ^m,n+l ; ^m,n+l ) ^i.'^mn^^nin-; ^mii) 

and similarly for t{x^ t, u). The quantities Vi of eq. I|126(l can be chosen to be 

Vl = Xm+l,n V2 = tm+l.n ^3 = Xm,n+1, /l Q0^ 

V4 — Em,n+l7 — Um,n+l- 

However, in (|131(l Um+i,n and Um.n+i cannot be expressed in terms of Umn, 
since eq. H114(l also involves Um-i,n- Differentiating (|131(l with respect to e.g. 
Wjnn- we find that ^ cannot depend on u: 

d^{Xmn , tmn , U„in ) ^ (133) 

Since we have i„+i„ — tmn and x,nn+i — the two equations (|131|l yield 

OXrn.n. C/lrn.n. 
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respectively. The same is obtained for the coefficient r, so finally we have 



? = eo, T = To, (135) 

where and tq are constants. 

Now let us apply pr X to eq. (|114|) . We obtain 

4'mn+l — 4>mn — Ti — ^{'Pm+ln — '2,(j)mn + 0m-l,n) — 0. (136) 

In more detail, eliminating the quantities Va in eq 1132|) we have 



+ 4'{Xmn - hi,tmn,Um~l,n)] = 0. (137) 



We differentiate eq. H137() twice, with respect to Um+i,n and itm-i,n respectively. 
We obtain 



- 0, (138) 
that is 

4'nin — ^(^mn, im^n)^ni7i ~t~ B(^Xm.n: ^m^n)- (139) 

We substitute of eq. (|139(l back into eq. 1)137(1 and set the coefficients 

of Um+i,m Umn, Um-in and 1 cqual to zero separately. From the resulting 
determining equations we find that A(xmn,tmn) = must be constant and 
that B{x,t) must satisfy the discrete heat equation H114I) . The result is that the 
symmetry algebra of the system (|114f) - (|116|l is very restricted. It is generated 

by 

Po^dt, Pi=d^, W = udu, S^S{x,t)du (140) 

and reflects only the linearity of the system and the fact that it is autonomous. 

The dilations, expansions and Galilei boosts, generated by D, K and B of 
eq. H129|l in the continuous case are absent on the lattice H115|l and H116|l . Other 
lattices will allow other symmetries. 

3.3.3 Discrete Heat Equation Invariant Under Dilations 

Let us now consider a five-point lattice that is also uniform and orthogonal. We 
put 

Um,n+1 ^m.n ^m+ln '^'^ra.n "t" ^m — In (141) 

^m+ln 2x77i7i -\- Xjn — lji — Xmn+1 ^mn — (142) 

tm+ln tmn — ^m,n+l ^^rn^n ^m.n— 1 — 0. (143) 
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The variables Va that we shall substitute from eq. p41ll . (|142|l and p43|l are 

Xm+i,n, Un+i,n, Xm,n+i, tm,n+i and UjriM+i- Applying pr X to eq. H142() we 
obtain 

^i^^mn^ '^^m^n ^mn—l: ^m,n+l) C(*^mm ^mn; ^mn) 0- (-^^^) 

In eq. H145() Umn and Um,n+i are independent. Differentiating with respect to 
Wmn we find dS^mn I dumn — and hence ^ does not depend on u. Differentiating 
p45|) with respect to tmn-i we obtain d£,m,7i+i/9tm,n+i — 0. Thus, ^ depends 
on X alone. Eq. (|f44|l can then be solved and we find that ^ is linear in x. 
Applying prX to eq. p43|) we obtain similar results for T{x,t,u). Finally, 
invariance of the lattice equations (|f 42|l and (|f 43|l implies: 

^ ^ ax + b, T = ct + d. (146) 

Let us now apply pr X to eq. (|f 4I|I . We obtain, after using the PAS ()I4I|I - 



^4^mn ~H ^m— In 



+ (2a-c) . ^2 =0- (147) 

Notice that Um,n+i (and hence 0m,n+i) depends on Um+i,n and Um-i,n, whereas 
all terms in eq. (|f 47|l depend on at most one of these quantities. Taking the 
second derivative dum+i,ndum-i,n of eq. (|147|l . we find 

= 0, i.e. (j) = A{x, t)u + B{x, t). (148) 

We substitute this expression back into (|147|l and find 

A{x, t) ^ Ao^ const (149) 

and see that B{x,t) must satisfy the system (|141ll - (|143|l . Moreover, we find 
c = 2a in eq. 1] 146(1 . Finally, the symmetry algebra has the basis 

Po = dt, Pi=a„ W^udu, D = xd^ + 2tdt (150) 
S = S{x,t)du. (151) 

Thus, dilational invariance is recovered, not however Galilei invariance. Other 
symmetries can be recovered on other lattices. 

3.4 Lorentz Invariant Difference Schemes 
3.4.1 The Continuous Case 

Let us consider the PDE 

Uxx - utt = 4/(w). (152) 
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Eq. H152|) is invariant under the Poincare group of 1 + 1 dimensional Minkowski 
space for any function f{u). Its Lie algebra is represented by 

Po = dt, Pi =d^, L = td^ + xdt. (153) 

For specific interactions f(u) the symmetry algebra may be larger, in particular 
for f = e"^, f — , or f = au + (3. 

Before presenting a discrete version of eq. (|152() . we find it convenient to 
pass over to light cone coordinates 

y^x + t, z = x-t (154) 

in which eq. 1152|l is rewritten as 



■yz 



f{u) (155) 



and the Poincare symmetry algebra (|153|l is 

Pi=dy, P2^dz, L^ydy-zdz. (156) 

3.4.2 A Discrete Lorentz Invariant Scheme 

A particular PAS that has eq. H155|) as its continuous limit is 

^m+ln+l '^m+ln ^mn \ 

w ^ T = f{Umn) (157) 

Vm+ln — 2y„in + Um-ln = ymn+1 — Umn = (158) 

To find the Lie point symmetries of this difference scheme, we put 

X = v{y^ u)dy + ^{y, z, u)dz + (l){y, z, u)du. (160) 

We apply the prolonged vector field pr x first to eq. p 58|) and (|I59II , eliminate 
ym+in, J/m,n+i, z^+in and z^.n+i, usiug the system p58|) . If59l) and notice that 
all Uik that figure in the obtained equations for -qik and are independent. 

The result that we obtain is that 77 and ^ must be independent of u and 
linear in y and z, respectively. Finally we obtain 

^ = aj/ + 7, r]^Pz + S (161) 

(a, . . . ,5 are constants). Invariance of eq. HI57() implies that the coefficient (p in 
the vector field HI60|I must be linear in u and moreover have the form 

(j)^ Au + B{y,z) (162) 

where A is a constant. Taking If 6f|) and If 62|) into account and applying pr A 
to eq. (|f57() . we obtain 

/ 4 r}\ r / \ , P'm-\-lTL-\-l Prnn+1 Pfn+ln ^ Pmn 
[A- a~ P)f(Umn) ^ 7 ^ w s 

+ Bmn)f'{Umn)- (163) 
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Differentiating eq. with respect to u,„„ we finally obtain the following 

determining equation: 

df (P f 

(a + P)-r^ + [Aumn + B{y„,n, z^^)]—^ = 0. (164) 

For f{umn) arbitrary, we find /3 = —a, A = B = Q. Thus for arbitrary f{u) 
the scheme (|157|) - (|158|) has the same symmetries as its continuous limit. The 
point symmetry algebra is given by eq. H15()(l . i.e. it generates, translations and 
Lorentz transformations. 

Now let us find special cases of f{u) when further symmetries exist. That 
means that eq. p64l) must be solved in a nontrivial manner. Let us restrict 
ourselves to the case when the interaction is nonlinear, i.e. 

'^^ /(m™,„)^0. (165) 



Then we must have 

B{ymn, Z,nn) ^ B ^ COUSt . (166) 

The equation to be solved for f{u) is actually eq. (|163|l which simplifies to 

{Au + B)f{u)^{A-a-f3)f{u). (167) 
For A^ the solution of eq. H167|) is 

/ = fou^ (168) 

and the symmetry is 

2 

Di = ydy + zdz udu. (169) 

P- 1 

For A = 0, i? 7^ we obtain 

/ = /oeP" (170) 

and the additional symmetry is 

D2 = ydy + zd, - 2du. (171) 

Thus, for nonhncar interactions f{u), f" ^ 0, the PAS (|157|l - H159|l has 
exactly the same point symmetries as its continuous limit (|155|) . 
The linear case 

f{u)^Ru + T (172) 

is different. The PDE (|155|l in this case is conformally invariant. This infinite 
dimensional symmetry algebra is not present for the discrete case considered in 
this section. 
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^ — o — o — o — o- 

Un-1 Un Un+1 

Figure 3: A monoatomic chain 

4 Symmetries of Discrete Dynamical Systems 
4.1 General Formalism 

In this chapter we shall discuss differential-difference equations on a fixed one- 
dimensional lattice. Thus, time t will be a continuous variable, n e Z a dis- 
crete one. We will be modeling discrete monoatomic or diatomic molecular 
chains with equally spaced rest positions. The individual atoms will be vibrat- 
ing around their rest positions. For monoatomic chains the actual position of 
the n-th atom is described by one continuous variable Un{t). For diatomic atoms 
there will be two such functions, u„(i) and Vn(t). 

Only nearest neighbour interaction will be considered. The interaction are 
described by a priori unspecified functions. Our aim is to classify these functions 
according to their symmetries. 

Three different models have been studied |42[l23ll5^ . They correspond to 
Fig. 13 21 and [3 respectively. 

The model illustrated on Fig. |2| corresponds to the equation 03 

u„(t) - i^„(i,u„_i(t),u„(t),u„+i(t)) =0. (173) 

Fig. 21 could correspond to a very primitive model of the DNA molecule. The 
equations are p51 

Un = Fn{t,Un~l{t),Un{t),Un+l{t),Vn-l{t),Vn{t),Vn+lit)) ^ 

Vn = Gn{t,U„-.i{t),U„{t),Un+lit),Vn-lit),Vnit),Vn+lit)) ^ 0. (174) 

The model corresponding to Fig. El already took translational and Galilei 
invariance into account, so the equations were 

Un = Fn i^n , t) + G„ (?/„-!, t) 

in = Kn{^n,t)+Pn{Vn~l,t) (175) 
^n — Vn -^n 1 — '^n+l Un- 

Dissipation was ignored in all three cases, so no first derivatives are present. 
In these lectures we shall only treat the case H173f) . The lattice is fixed, i.e. 
it is given by the relation 

Xn — hn -f xo (176) 

with h and xq given constants. With no loss of generality we can choose h — \, 
xo = 0, so that we have Xn — n. 
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Un-1 Un U„+l 
Vn-1 Vn Vn+i 



Figure 4: A diatomic molecule with two types of atoms on parallel chains 



Vn+1 Un+1 



Figure 5: A diatomic molecule with two types of atoms alternating along one 
chain 



Our aim is to find all functions Fn for which eq. (|173() allows a nontriv- 
ial group of local Lie point transformations. We shall also assume that the 
interaction is nonlinear and that it does indeed couple neighbouring states. 

Let us sum up the conditions imposed on the model H173|) and on the sym- 
metry studies. 

1. The lattice is fixed and regular (a;„ = n). 

2. The interaction involves nearest neighbours only, is nonlinear and coupled, 
d^Fn ,^ dFn ,^ dFn 



duiduk 



^0, 



dUn-l 



du 



^0. 



n+l 



(177) 



3. We consider point symmetries only. Since the lattice is fixed, the trans- 
formations are generated by vector fields of the form |38j 



X = T{t)dt + 4>nit, Un)dUn 



(178) 



We also assume that r(i) is an analytical hmction of t and (j)n{t,Un) is 
also analytic as a function of t and u„. 

The symmetry algorithm is the usual one, namely 

wXEn |b„=o= 0. (179) 

The prolongation in eq. (| 179(1 involves a prolongation to ^-derivatives and 
Um and to all values of n figuring in eq. I|173|l . i.e. n ± 1. 
The terms in the prolongation that we actually need are 



n+l 



(180) 



k—n—1 
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The coefficient is calculated using the formulas of Chaptern](or e.g. Ref. |48| ) . 
We have 

0** = D^<j)n - {DjT)Un - 2iDtT)il„. (181) 

Applying pr^^^ X to eq. (|173|l and replacing il from that equation, we get an 
expression involving (wn)^, {unY ^^'^ (un)"- The coefficients of all of 

these terms must vanish separately. The first three of these equations do not 
depend on Fn and can be solved easily. They imply 



(j>n{t, Ur, 



T{t) + ar 



T(t), d„ 



0. 



(182) 



The remaining determining equation is 



1 



-TtttUn + f3n + [ a-n + -ZT' ] F, 



-tF„ 



2 

E 



1 



and the vector fields realizing the symmetry algebra are 



X = T{t)dt 



1 X 

-T(t) +a. 



F„ 



dUn- 



(183) 



(184) 



Since we are classifying the interactions Fm we must decide which functions 
Fn will be considered to be equivalent. To do this we introduce a group of 
"allowed transformations", or a "classifying group". We define this to be a 
group of fiber preserving point transformations 

M„(i) = ri„(u„(t),t,5), i^i{t,g), n = n, (185) 

taking eq. (|173|l into an equation of the same form 

Un{i) = F„(t,t2„_i(t),{t„(t),{t„+i(i)) = 0. (186) 

That is, the allowed transformations can change the function F^ (as opposed 
to symmetry transformations), but cannot introduce first derivatives, nor other 
than nearest neighbour terms. These conditions narrow down the transforma- 
tions (|185|l to linear ones of the form 

Unit) = 4^u„(t) + B„(i), i^i{t) 
Vtt 

An,t = 0, it ^ 0, An 7^ 0, n = n. 



(187) 



Eq. (|173|l is transformed into 

<,tt = An^{tty^/'^\ Fn{t, Un-l,Un,Un+l) 



-^^Aniur^/'iiur 



An 



Un{i)~ Bn,tt>. (188) 
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The transformed vector field p84|l is 



X ^ T{t)tt{t)di + 



dUn. (189) 



In eq. H188|l and H189|l r(i), a„, /3„ and -F„ are given, whereas i{t), An and Bn{t) 
are ours to choose. We use these quantities to simphfy the vector field X. 

Our classification strategy is the following one. We first classify one-dimen- 
sional subalgebras. Thus, we have one vector field of the form (I184|l . If T{t) 
satisfies T(t) ^ in some open interval, we use t{t) to normalize T(t) — 1 and 
Bn{t) to transform /3„(i) into /3„(t) = 0. If we have r(t) = 0, a„ 7^ 0, we use 
Bn{t) to annul f3n{t)- The last possibility is r(i) = 0, a„ = 0, (3n{t) 7^ 0. Then 
we cannot simplify further. The same transformations will also simplify the 
determining equation 11831) and we can, in each case, solve it for the interaction 

Once all interactions allowing one dimensional symmetry algebras are de- 
termined, we proceed further structurally. We first find all Abclian symme- 
try groups and the corresponding interactions allowing them. We run through 
our list of onc-dimensional algebras and take them in an already establlished 
"canonical" form. Let us call this element Xi (in each case). We then find all 
elements X of the form H184(l that satisfy [Xi, X] = 0. We classify the obtained 
operators X under the action of a subgroup of the group of allowed transforma- 
tions, namely the isotropy group of Xi (the group that leaves the subalgebra 
Xi invariant). For each Abelian group we find the invariant interaction. 

^From Abelian symmetry algebras we proceed to nilpotent ones, then to 
solvable ones and finally to nonsolvable ones. These can be semisimple, or they 
may have a nontrivial Levi decomposition. 

All details can be found in the original article 021 > here we shall present the 
main results. 



4.2 One-Dimensional Symmetry Algebras 

Three classes of one-dimensional symmetry algebras exist. Together with their 
invariant interactions, they can be represented by 



^1, 



X = dt+ anUndUn 

Fn{t,Uk) = /„(C„-l,Cn,^«+l)e°"* 



Uke 



fc = n— l,n,ri + l. 



(190) 



X = a„u„5u„ 



(191) 
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Fn{t,Uk) = -T^Un + /n(i,Cn-l7?n+l) 



(192) 



We see that the existence of a one-dimensional Lie algebra implies that the 
interaction F is an arbitrary function of three variables, rather than the original 
four. The actual form of the interaction in eq. H190|l . (|191|l and (|192() was 
obtained by solving eq. 1183|l . once the canonical form of vector field X in 
eq. (|190|l . p91|l . or (|192(l was taken into account. 

4.3 Abelian Lie Algebras of Dimension N > 2 

Without proof we state several theorems. 

Theorem 1. An Abelian symmetry algebra of eq. (|173|l can have dimension N 
satisfying I < N < A. 

Comment: For = 1 these are the algebras A1.2 and Ai,^ of eq. H19U|) . 

(dnH) and fl^ . 

Theorem 2. Five distinct classes of interactions Fn exist having symmetry 
algebras of dimension N = 2. For four of them the interaction will involve an 
arbitrary function of two variables, for the fifth a function of three variables. 

The five classes can be represented by the following algebras and interactions. 

^2,1 : Xi^ dt+ ainUndUn, X2 = a2nUndUn 




k = 



nil 



(193) 



^2,2 : 



Xi ^ dt+ OnUndUn X2 = e°"'^dUn 
Fn = a^Un + e""*/«(Cn-l,Cn+l) 

6 = Ufee-"^*-u„e-""*, k^n±l 



(194) 



^2,3 : 




(195) 



^2,4 : 




n 



(196) 
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p _ {PlnP2n ~ /3ln/?2n)Ur!.+ l " {Pln+lP2n — /3ln/32n+l) 

i — (/3ln/32n+l — /3ln+l/32n)u„-l + (/3l„+l/32n-l — Pln-lP2n+l)Un 
+ (/3ln-l/32n — /3ln/32n-l)u„+l 

-Fn = Ck^Uk~Un, fc = n±l (197) 

The algebra ^2,5 is of particular physical significance since Xi and X2 in 
eq. H197(l correspond to translation and Galilei invariance for the considered 
chain. Unless we are considering a molecular chain in some external field, or 
unless some external geometry is imposed, the symmetry algebra ^2,5 should 
always be present, possibly as a subalgebra of a larger symmetry algebra. 

Theorem 3. Precisely four classes of three-dimensional symmetry algebras ex- 
ist. Only one of them contains the A2^5 subalgebra and can be presented as 

^3,4 Xi = dUn, X2 = tdUn, 

X3 ^ I3n{t)dun, /3„+i ^ /3„, ^„ ^ 0. (198) 
The invariant interaction is 

Fn = . ^" . {Un+l-Un) + fnit,0, (199) 

Pn+1 — Pn 

i = (/3„ - /3„+i)m„_i + (/3„+i - /3„-l)u„ + (/3n-l - /3n)Wn+l- (200) 

For ^3^1, A3. 2 and A3. 4 see the original article j42j . 

Theorem 4. There exist precisely two classes of interactions in eq. (|173|l 
satisfying conditions H177II . allowing four-dimensional symmetry algebras. Only 
one of them contains the subalgebra ^2.5 and is represented by the following. 

In - 7«+l 

Xi^dun, X2 = tdun, X3 = -01 (i)7„(9u„, (201) 

Xi = 4'2{t)jndUn 

7n+i7^7«, 7« = 0, -01 -02 - -01 ^/'2 = const 7^ 
with ^ as in eq. \2UU\) and tpi , ■02 satisfying 

B{t)i;,^Q, i = l,2 

4.4 Some Results on the Structure of Lie Algebras 

Let us recall some basic properties of finite-dimensional Lie algebras. Consider 
a Lie algebra L ~ {Xi,X2, . . . ,Xn}, where the elements Xi form a basis. To 
each algebra L one associates two series of subalgebras. 
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The derived series consist of the algebras 

L° = L, L^ = DL= [L, L], = D'^L = [DL, DL], . . . 

= = (202) 

The algebra of commutators DL is called the derived algebra. If we have DL = 
L, the algebra L is called perfect. If an integer A'' exists for which we have 
D^L = {0}, the algebra L is called solvable. 
The central series consist of the algebras 

Lo = L, Li=L^ = [L,L], L2 = [L,Li],...Ln=[L,Ln-i],... (203) 

If there exists an integer N for which we have L^ = {0}, the algebra L is 
called nilpotent. Clearly, every nilpotent algebra is solvable, but the converse is 
not true. 

Let us consider two examples 

1. The Lie algebra of the Euclidean group of a plane: e(2) ~ {L3,Pi,P2}. 
The commmutation relations are 

[L3,Pi]=P2, [L3,P2]=-Pi, [Fi,P2]= 0. (204) 

The derived series is 

L = {L3,Pi,P2}D DL = {Pi,P2}, D^L = {0} 

and the central series is 

LdLi = {Pi,P2} = L2 = L3 = ... 

Hence e(2) is solvable but not nilpotent. 

2. The Heisenberg algebra Hi ~ {Xi,X2,X3} where the basis can be re- 
allized e.g. by the derivative operator, the coordinate x and the identity 
1: 

Xi=d^,X2=X,X3 = 1. 

We have 

[Xi, X2] = X3, [Xi, X3] = [X2, X3] = (205) 

and hence 

DL = {X3}, D^L = 0. 
Li = {X3}, L2 = 0. 

We see that the Heisenberg algebra is nilpotent (and solvable). 
An Abelian Lie algebra is of course also nilpotent. 

We shall need some results concerning nilpotent Lie algebras (by nilpotent 
we mean nilpotent non- Abelian) . 
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1. Nilpotent Lie algebras always contain Abelian ideals. 

2. All nilpotent Lie algebras contain the three-dimensional Heisenberg alge- 
bra as a subalgebra. 

We shall also use some basic properties of solvable Lie algebras, where by 
solvable we mean solvable, non- nilpotent. 

1. Every solvable Lie algebra L contains a unique maximal nilpotent ideal 
called the nilradical NR{L). The dimension of the nilradical satisfies 

i dim(L) < dim7Vi?(L) < dim(L) - 1. (206) 

2. If the nilradical NR{L) is Abelian, then we can choose a basis for L in 
the form {A"i, . . . , Xn, Yi, . . . , i^n}, m < n, with commutation relations 

[X„Xk]=0, [X,,Yk] = {Ak)^JX,, [Y,,Yk] ^ ClkXi. (207) 

The matrices Ak commute and are linearly nilindependent (i.e. no linear 
combination of them is a nilpotent matrix). 

If a Lie algebra L is not solvable, it can be simple, semisimple, or it may have 
a nontrivial Levi decomposition |32| . A simple Lie algebra L has no nontrivial 
ideals, i.e. 

IQL, [1,1] Q I, [L,I]CI (208) 

implies / ^ {0}, or / = L. 

A semisimple Lie algebra L is a direct sum of simple Lie algebras Li 

Lr^ Li®L2® ...®Lp, [L^,Lk]^ 0. (209) 

If L is not simple, semisimple, or solvable, then it allows a unique Levi decom- 
position into a semidirect sum 

[5, 5] =5, [R,R]ciR, [S,R]CR (210) 

where S is semisimple and R is solvable; R is called the radical of L, i.e. the 
maximal solvable ideal. 

Let us now return to the symmetry classification of discrete dynamical sys- 
tems. 

4.5 Nilpotent Non-Abelian Symmetry Algebras 

Since every nilpotent Lie algebra contains the three-dimensional Heisenberg al- 
gebra, we start by constructing this algebra. Hi ~ {Xi, ATa}. The central 
element X^ of eq. ()205|l is uniquely defined. We start from this element, take it 
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in one of the standard forms p9U|) . (|191|l . or p92|) . then construct the two com- 
plementary elements Xi and X2 ■ The result is that two inequivalent realizations 
of Hi, exist namely: 

Fn ^ fniS,n+l,Cn-l), S,k = Uk - Un, k = n±l (211) 



Xz^ {t + 7„)e''"*au„, d„ = 0, 7„ = 0, 7„+i ^ 7„ 
7„+i - 7„ 



7n+l - 7n 

C = (7™ - 7„+i)w„_ie"""-i* + (7„+i - 7„_i), Une"""* 

+ (7n-l - 7n)'^n+ie"''"+'*. 

Notice that N^ i contains the physically important subalgebra ^2^5. whereas 
■^3,2 does not. 

Extending the algebras N^^i and by further elements, we find that N^^i 
gives rise to two five-dimensional nilpotent symmetry algebras Nc,^^. and iV3^2 to 
a four-dimensional one N^^i. 

Here we shall only give A^5,i and A^5,2 which contain N^^i and hence ^2,5: 

^5,fc : Xi = dun, X2 = tdun, X3 = i h 7„ j dun, 

^4 = (^^^^^ + rnt^ dUn, X^^dt, k^l,2 (213) 

2(k — 1) 

Fn - ^ -iUn+1 - Un) + fniO 

7n+l - In 

with ^ as in eq. (|212|l . 

4.6 Solvable Symmetry Algebras with Non-Abelian Nil- 
radicals 

We already know all nilpotent symmetry algebras, so we can start from the 
nilradical and extend it by further non-nilpotent elements. The result can be 
stated as a Theeorem. 

Theorem 5. Seven classes of solvable symmetry algebras with non-Abelian nil- 
radicals exist for eq. (|173|) . Four of them have N3 1 as nilradical, three have 

For N^^i we can add just one further element Y , namely one of the following 
SN4,i : Y = tdt + (l- + a]undun, a 7^ — 
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F„ ^ K+i - u„)e(''-^/^)/('^+i/^)/„(0 (214) 

SN4^2 ■■ Y = tdt + (2u„ + t'^)dun 

K =ln(w„+i-M„)+/„(0 (215) 

SN4^3 : Y = Undun 

Fn = {Un+l-Un)fn{0- (216) 

In all above cases we have 

u„^i-7/„ ^ (217) 

SNi^i : Y = tdt+ Jndun, 7,1+1 7^ 7„, 7„ = 

J^.=expf-2 ""+^'"" V»(0 (218) 



7n+l - 7n , 
C = (7n - 7n+l)Wn-l + (in+l " 7«-l)"n 

+(7„_i - 7„)m„+i. (219) 
For iVs^i we can also add at most one non-nilpotent element and we obtain 

SNe^i : Y^tdt+(^^ + a^ u^du„ 

F„=C„C("-3/2)/(a+l/2)^ ^^_1^ ^_,3 ^220) 

SNe^2: Y = tdt + [2u.n + (a + b-fr,)t^]dur„ + ^ 

=c„ + (a + 67„)ln^ (221) 

SNq^S: Y^tdt+PndUn, Pn^^A + Bjn, 7^ 

Fn = c„e'^ (222) 
. 

(7n ~ ln+l)pn-l + {Sn+1 - ln-l)Pn + {in-l - ln)pn+l 

In all cases ^ is as in eq. (j219|l . 

4.7 Solvable Symmetry Algebras with Abelian Nilradicals 

The results in this case are very rich. There exist 31 such symmetry algebras 
and their dimensions satisfy 2 < d < 5. 

For all details and a full complete list of results we refer to the original 
article. Here we give just one example of a five-dimensional Lie algebra with 
NR{L) = ^4,1. 

5*^54 : Xi = dun, X2 = tdun, X3 = e*7„9Mn, 
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e *7„9u„ 

Y^dt + aundun a^O, 7„^7„+i, 7„ = 

F„ = ^"(^"+1-"") + e'^VnlO (223) 

7n+l - 7n 

C = [(7n - 7n+l)'"n-l + iln+1 - 7n-l)"« 
+ (7ri-l - 7ri)u,i+l]e"'''. 

4.8 Nonsolvable Symmetry Algebras 

A Lie algebra that is not solvable must have simple subalgebra. The only simple 
algebra that can be constructed out of vector fields of the form (|184|l is sl(2, R). 
The algebra and the corresponding invariant interaction can be represented as: 

5*3,1 : Xi = dt, X2 = tdt + ^u„i9u« 

A3 =t^dt + tundun (224) 

-^n o"/n (Cn— 1 ; Cn+l); — • 

This algebra can be further extended to a four, five or seven-dimensional sym- 
metry algebra. In two cases the algebra will have an ^2,5 subalgebra, namely 
A'S'5,i: In addition to Ai, A2, A3 of we have 

A4 = 9u„, A5 = tduj^ 

F„ = (u„+i-u„)-V„(0, ^^^hl±l^, (225) 

^n — 1 

NSt.i: The additional elements are 

A„ = dUn, A5 = tdUn, Xe = JndUn, A7 = t^nOUn 

7„+i ^ 7„, 7„ = 0. 
The invariant interaction is 



(226) 



Fn = s„[(7„ - 7„+i)w„_i + (7„+i - 7„_i)u„ 

+ (7n-l -7n)"n+l]"^ S„ = 0, S„ 7^ 0. (227) 

4.9 Final Comments on the Classification 

Let us first of all sum up the discrete dynamical systems of the type p73|l with 
the largest symmetry algebras 
We put 

C = (7« - 7n+l)Wn-l + (7n+l " 7n-l)"n + iln-1 " 7n)Wn+l (228) 

and find this variable is involved in all cases with 7, or 6-dimensional symmetry 
algebras. 
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The algebras and interactions are given in eq. H226|) , (|220l) , (|221|l and H222|) , 
respectively. 

A natural question to ask is: Where is the Toda lattice in this classification? 
The Toda lattice is described by the equation 

u„,« = e""-^-"" - e""-""+^ (229) 

This equation is of the form (|173|l . It is integrable [S7] and has many interesting 
properties. In our classification it appears as a special case of the algebra 
i.e. 

ii„ = expf-2 ""+^~^" Vn(0, (230) 

with 

/n(C) = -l + e«/^ 7« = 2n. (231) 

Thus, its symmetry group is four-dimensional. We see that the Toda lattice is 
not particularly distinguished by its point symmetries: other interactions have 
larger symmetry groups. Even in the SN^^^ class two functions have to be 
speciahzed (see eq. ()231ll to reduce H230|l to (|229|l . 



5 Generalized Point Symmetries of Linear and 
Linearizable Systems 

5.1 Umbral Calculus 

In this chapter we take a different point of view than in the previous ones. In- 
stead of purely point symmetries, we shall consider a specific class of generalized 
symmetries of difference equations that we shall call "generalised point symme- 
tries" . They act simultaneously at several, or even infinitely many points of a 
lattice, but they reduce to point symmetries of a differential equation in the 
continuous limit. 

The approach that we shall discuss here is at this stage applicable either to 
linear difference equations, or to nonlinear equations that can be linearized by 
a transformation of variables (not necessarily only point transformations). 

The mathematical basis for this type of study is the so called "umbral cal- 
culus" reviewed in recent books and articles by G.G. Rota and his collabora- 
tors |55I54I1(J] . Umbral calculus provides a unified basis for studying symmetries 
of linear differential and difference equations. 

First of all, let us introduce several fundamental concepts. 

Definition 1. A shift operator Ts is a linear operator acting on polynomials or 
formal power series in the following manner 

Tsfix) ^ f{x + 6), xeR, SeK. (232) 
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For functions of several variables we introduce shift operators in the same 
manner 

= f{xi,...,Xi-i,Xi+Si,Xi+i,...,Xn). (233) 

In this section we restrict the exposition to the case of one real variable 
X G R. The extension to n variables and other fields is obvious. 

Definition 2. An operator U is called a "delta operator" if it satisfies the 
following properties 

1) It is shift invariant; 

TsU = UTs, V5 e R (234) 

2) 

Ux^^c^O, c = const (235) 

3) 

Ua = 0, Va (236) 
and the kernal of U consists precisely of all constant. 
Important properties of delta operator are: 

1. For every delta operator U there exists a unique series of basic polynomials 
{Pnix)} satisfying 

Pq{x) = 1 p„(0) =0, n > 1, Upn{x) = np„_i(x). (237) 

2. For every umbral operator U there exists a conjugate operator /3, such 
that 

[U,xl3] = \. (238) 

The operator /3 satisfies 

(3 = {U)-^, U=[U,x\. (239) 

The expression 

U=U *x=[U, x] (240) 

is called the "Pincherle derivative" of U [CTI^ITIl] . 

For us the fundamental fact is that the pair of operators, U and x/?, satisfies 
the Heisenberg relation H238|) . 

Before going further, let us give the two simplest possible examples. 

Example 6. The (continuous) derivative 

^ = 5., /3=1 

P0 = 1, Pl=X, ...,Pn=x",... ^ ' 
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Example 7. The right discrete derivative 

U = A+ = ^, (3 = T-^ 
Po = l, Pi=x, P2 = x{x^S) (242) 
P„ ^ x{x - S) . . . {x - {n - 1)S). 

For any operator U one can construct /3 and the basic series will be 

P„ = (x/?)" -1, n e N. (243) 

5.2 Umbral Calculus and Linear Difference Equations 

First of all, let us consider a Lie algebra L realized by vector fields 

Xa = fa{xi,...,X„)dXa (244) 

[Xa,X,]^C2kX,. (245) 

The Heisenberg relation H238|l allows us to realize the same abstract Lie algebra 
by difference operators 

X^ ^ faixif3i,X2f32,...,XnPn)A^^, [A^^ , XaPa] = 1 , a=l,...n. (246) 

As long as the functions fa are polynomials, or at least formal power series in 
the variables Xa, the substitution 

Xa XaPa, A^^ (247) 

preserves the commutation relations H245(l . 

We shall call the substitution (|247|l and more generally any substitution 

{U,p}^{Uj} (248) 

an "umbral correspondence" . This correspondence will also take the set of basic 
polynomials related to {U, /3} into the set related to the pair {u,(3}. 

We shall use two types of delta operators. The first is simply the derivative 
U — dx, for which we have /3 = 1. The second is a general difference operator 
U — A that has dx is its continuous limit. We put 

A = - ^ OfeTl /, m G Z, l<m (249) 



k=l 



where ak and 5 are real constants and Ts is a shift operator as in eq. 1)232(1 . 
Condition (|234ll is satisfied. Condition 123()l) implies 

m 

^Ofc^O. (250) 

k=l 
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We also require that for 5 ^ 0, we should have A ^ dx- This requires a further 
restriction on the coefficients ak, namely 



^afc/s = l. (251) 



k=l 



Then relation (|235|l is also satisfied, with c = 1. 
More generally, we have, for A as in H249|l 



k=l 



We define 



and thus 



7g 



q=0 ^' k=l 



^^Gkk'' qeZ (252) 

k=l 



9=2 

Thus A goes to the derivative at least to the order 5. We can also impose 
7^ = 0, g = 2,3,...m-L (254) 

Then we have 

A = -^ + 0((5"-i). 

dx 

Definition 3. A difference operator of degree m — Hs a delta operator of the 
form 

^ m 

U = A = -J2akT^, (255) 

k=l 

where ak and 6 are constants, Tg is a shift operator and we have 

m ni 

J2ak = 0, ^flfcfc-l. (256) 

k=l k=l 

Comment: A — A is a difference operator of the same degree as A. 
Theorem 6. The operator (3 conjugate to A — (1/(5) X^fcL; o,kTg is 

/3=(^afcA;TM . (257) 
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Proof. /3= (A')-^ = [A,a;]-i 

^ / m m 

[A,x] = -(y.ak{x + k5)T^ -xY^auT 

^k=l k=l 



□ 



Examples: 



A3 = -_(r2-6T + 3 + 2r-i), 13 =i ^ j (259) 

Comment: 

A^ = ^ + 0{S') A3 = A + 0(<53). 
ox ox 

Now let us apply the above considerations to the study of symmetries of 
linear difference equations. 

Definition 4. An umbral equation of order n is an operator equation of the 
form 

n 

J2&k{xP)A''f = g (260) 

fc=0 

where ak{x(3) and g{x(3) are given formal power series in xP and f{x(3) is the 
unknown operator function. 

For IS. — P = I this is a differential equation. For A as in (|249|l eq. (|260|l 
is an operator equation. Applying both sides of eq. H2t)U|l to 1 we get a difference 
equation. Its solution is 

f{x) = fixP) ■ 1. (261) 
More generally, an umbral equation in P variables is 

ni . ..Up 

^ aki...kpixi(3i,x2l32, ■ ■ .,XpPp)Agl,. . . A^^/(xi/3i, . . .,XpPp) 

p 

= g{xipi . . .XpPp), ^ni = n. (262) 

Example of an umbral equation 

A/ = a/, a^O. (263) 
(i) Take A^d^^ f{x) = e"'-' 
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(ii) Take A = A+ = P = T"! 

f{x + d)-f{x)^aSf{x). 

Take f{x) = A^: 

= aSX'' ^ A = (1 + aSy/^. 
We get a single "umbral" solution 

hix)^{i + asy-/'. 

The umbral correspondence gives: 

/2(x) = e^^^'" ■ 1. 



(264) 



(265) 



(266) 



If we expand into power series, we obtain fi{x) — f2{x), and of course we 
have 

lim/i,2(x) =e"^ 

6 — >0 

(iii) For comparison, take A ^ A" = (T - T-^)/25, (3^[{T + T-^)/2]-'^ 

fix + S)- f{x-S)^2Saf{x). (267) 
Putting f{x) ~ we get two values of A and 



/ = Ai{aS+ y/a^S^ + If^^ + A2 {a6 ~ y/a^P + lf/^ 

= Ai/1+^2/2. (268) 



We have 



lim/i(a;) = e"^ 



(269) 



but the limit of /2(a;) does not exist. The umbral correspondence yields 

'T + T-' 



fu{x) = exp 



ax 



1. 



Expanding into power series, we find fu — fi- The solution /2 is a non- 
umbral one. 

Theorem 7. Let A be a difference operator of order p. Then the linear um- 
bral equation of order n ]260(l has np linearly independent solutions, n of them 
umbral ones. 

There may be convergence problems for the formal series. 
Consider the exponential 



f{x) = e"^^ 



k=l 



(270) 
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For m~l > 3, /3 will involve infinitely many shifts i.e., each term in the expansion 
1)2 70(1 could involve infinitely many shifts. However 

P„(x) = {x(3r ■ 1 (271) 

is a well defined polynomial. For a proof see '36'. 

Let us assume that we know the solution of an umbral equation for A = 
and it has the form 

/(^) - E ^^-4^^"- (272) 
•■^ — ' n! 

ri=0 

Then for any difference operator A there will exist a corresponding umbral 
solution 

n=0 

where Pn{x) = {x(3)^ ■ 1 are the basic polynomials corresponding to A. 

5.3 Symmetries of Linear Umbral Equations 

Let us consider a linear differential equation 

Lu = 0, L= E ''k^,...,k,ixl,...,Xp)—^...-^. (274) 

The Lie point symmetries of eq. H213|) can be realized by evolutionary vector 
fields of the form 

p 

Q = (/."E^'".-.- (275) 
1=1 

The following theorem holds for these symmetries. 

Theorem 8. All Lie point symmetries for an ODE of order n > 3, or a PDE 
of order n > 2 are generated by evolutionary vector fields of the form H275|) with 
the characteristic Q satisfying 

Q^Xu + x{xi,...,Xp), (276) 

where x a solution of eq. \274\ and X is a linear operator 

p 

X = Y,U^i,---,Xp)dx^ (277) 

i=l 

satisfying 

[L,X] = \{x^,...,Xp)L, (278) 

i.e. commuting with L on the solution set of L. In eq. \27l^ A is an arbitrary 
function. 
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For a proof we refer to the hterature ^ . 

In other words, if the conditions of Theorem |S1 apply, then aU symmetries 
of eq. H213|) beyond those representing the hnear superposition principle, are 
generated by linear operators of the form (|277|l . commuting with L on the 
solution set of eq. (|274|l . 

Now let us turn to the umbral equation H2t)2|) with g = 0, i.e. 

J2 ak,...k,{xif3i,...Xpf3p)A'll...A';iu{xiPu...,Xp(3p)^ 0. (279) 

We shall realize the symmetries of eq. H279|) by evolutionary vector fields of 
the form 

p 

VE = QdOu, Qd ^ <t>D - iosA^u (280) 

i=l 

where 4>d and S,D,i are functions of Xil3i and u. The prolongation of ti^; will also 
act on the discrete derivatives A^*m. We are now considering transformations on 
a fixed (nontransforming) lattice. In the evolutionary formalism the transformed 
variables satisfy 

Xkfik = XkPk, Pk = Pk 
u{xkPk) = u{xkl3k) + XQD, |A|<1 (281) 

and we request that V' be a solution whenever ijj is one. 

The transformation of the discrete derivatives is given by 

Ai^i^u = A^^^^u + XA^^^^Q (282) 

etc., where A^^. are discrete total derivatives. 
In terms of the vector field H290|) we have 

= ^^QD, QT" = A.AfcQc (283) 

(we have put A^. = A^;. = A^). 

As in the continuous continuous case, we obtain determining equations by 
requiring 

V^VE{LDU)\Lni=^ (284) 

where Lou is the left hand side of eq. H279|l . 

The determining equations will be an umbral version of the determining 
equations in the continuous case, i.e. are obtained by the umbral correspondence 

^ Aj, Xi > Xif3i. 

The symmetries of the umbral equation H279|l will hence have the form (|280|l 
with 

Qd = Xdu + xixiPi, . . .XpPp) (285) 
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where Xd is a difference operator commuting witli Lu on the solutions of 
eq. (|279|l . Moreover, Xo is obtained from X by the umbral correspondence. 

We shall call such symmetries "generalized point symmetries" . Because of 
the presence of the operators (3i they are not really point symmetries. In the 
continuous limit they become point symmetries. 

Let us now consider some examples. 

5.4 Example of the Discrete Heat Equation 

The (continuous) linear heat equation in (1 + 1) dimensions is 

ut - u^x = 0. (286) 

Its symmetry group is of course well-known. Factoring out the infinite dimen- 
sional pseudo-group corresponding to the linear superposition principle we have 
a 6 dimensional symmetry group. We write its Lie algebra in evolutionary form 
as 



Po = utdu, Pi=Uxdu, W^udu 
B 

1 



K 



t^Ut + txux H — [x^ + 2t)u 
4 



du 



(287) 



where Pq, Pl^ B, D, K and W generate time and space translations, Galilei 
boosts, dilations "expansions" and the multiplication of w by a constant, re- 
spectively. 

A very natural discretization of eq. H286() is the discrete heat equation 



Atu- (Axfu^O, 



(288) 



where At and A^^ are each one of the difference operators considered in Sec- 
tion l5.2l We use the corresponding conjugate operators f3t and Px, respectively. 
The umbral correspondence gives us the symmetry algebra of eq. H288fl , starting 
from the algebra H287(l . Namely, we have 



pD 

3-D 



iAtu)du, Pf = {AxU)du, W = Udu 



= [2{tPt)AxU+(xl3x)u]du 



2tf3tAtu + xf3xAxU+ -u 



du 



{tptfAtU + {tPt){xl3x)AxU +-{ixPxf + 2tpt)u 



du. 



(289) 



In particular we can choose both At and A^; to be right derivatives 



At 



Tt-1 



l3t=T-\ A, 



- 1 



(290) 
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The characteristic Qk oi the clement then is 

Qk = Xku, Xk = (t' - <5tt)rr'At + txT-%-^A 



+ ^[ix'-Sx)T-' + 2tT,-\ (291) 

so it is not a point transformation: it involves u evaluated at several points. 
Each of the basis elements (|289|) (or any linear combination of them) provides 
a flow commuting with eq. 



ux = Xu. (292) 

Equations H288|l and (|292|l can be solved simultaneously and this will provide a 
difference analog of the separation of variables in PDEs and a tool for studying 
new types of special functions. 

5.5 The Discrete Burgers Equation and its Symmetries 
5.5.1 The Continuous Case 

The Burgers equation 

Ut = Uxx + 2mM:e (293) 

is the simplest equation that combines nonlinearity and dissipative effects. It 
is also the prototype of an equation linearizable by a coordinate transformation 
C-linearizable in Calogero's terminology |59| . 

We put u = Vx and obtain the potential Burgers equation for v: 

vt=vxx + vl. (294) 

Putting w — e" we find 

Wt = Wxx- (295) 

In other words, the usual Burgers equation (|293|l is linearized (into the heat 
equation (|295|l by the Colc-Hopf transformation 

u=^ (296) 

w 

(which is not a point transformation) . 

One possible way of viewing the Cole-Hopf transformation is that it provides 
a Lax pair for the Burgers equation: 

Wt = Wxx, Wx = UW. (297) 

Putting 

Wt — Aw, Wx — Bw, A — Ux + u^ , B = u 
we obtain the Burgers equation as a compatibility condition 

Ax-Bt + [A, B] = 0. (298) 
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Our aim is to discretize the Burgers equation in such a way as to preserve 
its hnearizability and also its five-dimensional Lie point symmetry algebra. We 
already know the symmetries of the discrete heat equation and we will use them 
to obtain the symmetry algebra of the discrete Burgers equation. This will be 
an indirect application of umbral calculus to a nonlinear equation. 



5.5.2 The Discrete Burgers Equation as a Compatibility Condition 

Let us write a discrete version of the pair H297() in the form: 

At(/. = A^^cj), A^(j) = u(j) (299) 

where we take 

A, = ^, ^- = ^V^- (300) 

Ot Ox 

The pair (|299|l can be rewritten as 

Atc/} ^ {AxU + uTxu)4}, Axcf) = u4>. (301) 

We have used the Leibnitz rule appropriate for the discrete derivative A^^ of 
H3UU|I . namely 

Ax/.g = fix)A,g + (r,.9)A,/. (302) 

Compatibility of eq. H301|l . i.e. AxAt4> = AtA^cj) yields the discrete Burgers 
equation 

Atu = , ^""^ tA:^(A^u + uT^u). (303) 

1 + dt[Aj:AxU + uT^u] 

In the continuous limit At d/dt, A^ d/dx, 1, (5^: = 0, (5t = we 

reobtain the Burgers equation H293|l |3()l28j . This is not a "naive" discretization 
like 

Atu = {A^fu + 2uA,^u (304) 
which would loose all integrability properties. 



5.5.3 Symmetries of the discrete Burgers Equation 

We are looking for "generalized point symmetries" on a fixed lattice. We write 
them in evolutionary form 

Xe = g(x, t, T^tS, T^TfA.u, T^T/Atu, . . .)du (305) 

and each symmetry will provide a commuting flow 

MA = Q- 

We shall use the Cole-Hopf transformation to transform the symmetry algebra 
of the discrete heat equation into that of the discrete Burgers equation. 
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All the symmetries of the discrete heat equation given in eq. (I289|l can be 
written as 

S^S{x,t,(j),T,,T,...) (306) 

where S' is a hnear operator (the same is true for any Unear difference equation) . 
For the discrete heat equation 

At0-(A,)20 = O (307) 

with At and A^; as in eq. H300|l we rewrite the flows corresponding to eq. (|289|l 
as 



2trr'A, + a;T-i + -4T-i 



2tr-iAt + xr-iA, 



1 



16 



X X 



(308) 



!^A6 = 



Let us first prove a general result. 

Theorem 9. Let eq. H306|l represent a symmetry of the discrete heat equation 
H3U7|) . Then the same linear operator S provides a symmetry of the discrete 
Burgers equation l|3U3l) . the flow of which is given by 



' S(f) 

WA = (1 + (5a;U)Aa; I — 



(309) 



where [S(j))/4i can he expressed in terms ofu{x,t). 

Proof. We request that eq. H3Q6(I and the Cole-Hopf transformation in eq. (|299|l 
be compatible 

^(A,0) = A,</.A. (310) 



I From here we obtain 



Ax^Scj)) — uS(j) 



(311) 



A direct calculation yields 



A.(^) =^[A,(50)-^(50)] 



(312) 



and H309|l follows. It is still necessary to show that S(p/4> depends only on u{x, t). 
The expressions for S(j) can be read off from eq. H308(l . From there we see that all 
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expressions involved can be expressed in terms of u{x,t) and its shifted values 
(using the Cole-Hopf transformation again). Indeed, we have 



Ax(j) = u(f), At(j) = w0 
T^cj) = {l + d:cu)(t>, Tf(j) = (1 + Stv] 
1 



^ 1 + SxU 



t: 



1 



where we define 



V = A^u + uTxU. 



(313) 



(314) 



Explicitly, eq. (|309|l maps the 6 dimensional symmetry algebra of the discrete 
heat equation into the 5 dimensional Lie algebra of the discrete Burgers equa- 
tion. The corresponding flows are 



WAi = {l + 5tv)Atu 

uxs = (l + 4w)A^ 

MA4 = (l+4w)Aj 

MA5 = (l + 4w)Ai 



2tTf 



V 

l + Stv 
1 



1 + SxU 



1 



1 + 4u 2 ^ l + 5xu 



-T7 



1 + Jtw ' 1 + 5tv 



1 



(315) 



1 + OxU 1 + OtV 



-T7 



1 .^-1. 1 



1 + 5xU ^ 1 + 5xU 
1 



l + Stv * 1 + (5tw 



1 



1 



2 ^ Vl + '^^w"' l + ^tw 



1 



WAe 



0. 



□ 



5.5.4 Symmetry Reduction for the Discrete Burgers Equation 

Symmetry reduction for continuous Burger equation: we add a compatible equa- 
tion to the Burgers equation 

Ut = Uxx + "iuUx 

u\ = Q{x,t,u,Ux^t) ^ (316) 
and solve the two equations simultaneously. Example: time translations. 

ux^ut^ 0. (317) 
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Then u — u{x) and 

Uxx + '2'UUx ^ Q ^ Ux + = K . (318) 
/,From here, we obtain three types of solutions 

u——, u = fc arctanhfcx, u = fc arctanfcx. (319) 

X 

Symmetry reduction in the discrete case. All flows have the form H309|l . The 
condition u\ = Q hence implies 

S(j} = K{t)(j), (320) 

where K{t) is an arbitrary function. This equation must be solved together 
with the discrete Burgers equation in order to obtain group invariant solutions. 

Let us here consider just one example, namely that of time translations, the 
first equation in (|315|l . Eq. (|320|l reduces to 

At^ = K{t)cl), (321) 

i.e. 

V ^ AxU + uTxU = K{t). (322) 
We rewrite the Burgers equation as 

Atu = — ^ AxV V = AxU + uTxU. (323) 
1 + Otv 

However, from H322|l we have v = K{t) and hence Atu — 0, K(t) — Kg — const. 
Since (j> satisfies the heat equation, we can rewrite (|321ll as 

Axxcl) = Ko(t). (324) 

The general solution of H324() is obtained for Ko ^ by putting cj) — and 
solving (|324|l for a. We find 

= ci(l + v^4)^/'^ + C2(l - ^/K'oSxr^'^■ (325) 
For Ko = we have 

(j) — ci + C2X. (326) 

Solutions of the discrete Burgers equation are obtained via the Cole-Hopf trans- 
formation 

u = — ■ — . (327) 

9 

The same procedure can be followed for all other symmetries. We obtain linear 
second order difference equations for (j) involving one variable only. However, 
the equations have variable coefficients and arc hard to solve. They can be reex- 
pressed as equations for u{x,t), again involving only one independent variable. 
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Thus, a reduction takes place, but it is not easy to solve the reduced equations 
explicitly. 

For instance, Galilei invariant solutions of the discrete Burgers equation must 
satisfy the ordinary difference equation 

2tT^u + x- K{t) + 2t6,,:uT,,:U + St (^IT^u + IS^uT^u 

+xuTxU + xAxU - + l^a; - K{t)[6xU + 6t{TxAxU i^"^^) 
+uT^u - uT^u) + T^uTlu + ^xWr^uT^w] = 
where t figures as a parameter. 
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